Transcriptome-wide Association Study Research Articles

Integrating meta‐analysis of genome‐wide association study with Pig Genotype‐Tissue Expression resources uncovers the genetic architecture for age at first farrowing in pigs

AbstractAge at first farrowing (AFF) is a reproductive trait with low heritability and high importance in the pig industry. To enhance the statistical power of genome‐wide association study (GWAS) and further explore the genetic nature of AFF, we first conducted GWAS meta‐analysis using three Yorkshire populations, and then integrated the Pig Genotype‐Tissue Expression (PigGTEx) resources to interpret their potential regulatory mechanism. Additionally, we compared the AFF in pig with the age at first birth (AFB) of human using GWAS summary statistics. We detected 18 independent variants in GWAS meta‐analysis and 8 genes in gene‐based association analysis significantly associated with AFF. By integrating the PigGTEx resource, we conducted transcriptome‐wide association study (TWAS) and colocalization analysis on 34 pig tissues. In TWAS, we detected 18 significant gene‐tissue pairs, such as DCAF6 in uterus and CREG1 in blood. In colocalization, we found 111 potential candidate tissue‐gene pairs, such as GJD4 and LYPLAL1. We found that the homologous gene, CHST10, might be the potential candidate gene between humans in AFB and pigs in AFF. In conclusion, integrating GWAS meta‐analysis and PigGTEx resources is a meaningful way to decipher the genetic architecture of complex traits. We found that DCAF6, CREG1, GJD4, and LYPLAL1 are candidate genes with high reliability for AFF in swine. The comparative analysis showed that CHST10 might play a potentially critical role in AFB/AFF across human and pigs.

Transcriptome-wide association mapping provides insights into the genetic basis and candidate genes governing flowering, maturity and seed weight in rice bean (Vigna umbellata)

BackgroundRice bean (Vigna umbellata), an underrated legume, adapts to diverse climatic conditions with the potential to support food and nutritional security worldwide. It is used as a vegetable, minor food crop and a fodder crop, being a rich source of proteins, minerals, and essential fatty acids. However, little effort has been made to decipher the genetic and molecular basis of various useful traits in this crop. Therefore, we considered three economically important traits i.e., flowering, maturity and seed weight of rice bean and identified the associated candidate genes employing an associative transcriptomics approach on 100 diverse genotypes out of 1800 evaluated rice bean accessions from the Indian National Genebank.ResultsThe transcriptomics-based genotyping of one-hundred diverse rice bean cultivars followed by pre-processing of genotypic data resulted in 49,271 filtered markers. The STRUCTURE, PCA and Neighbor-Joining clustering of 100 genotypes revealed three putative sub-populations. The marker-trait association analysis involving various genome-wide association study (GWAS) models revealed significant association of 82 markers on 48 transcripts for flowering, 26 markers on 22 transcripts for maturity and 22 markers on 21 transcripts for seed weight. The transcript annotation provided information on the putative candidate genes for the considered traits. The candidate genes identified for flowering include HSC80, P-II PsbX, phospholipid-transporting-ATPase-9, pectin-acetylesterase-8 and E3-ubiquitin-protein-ligase-RHG1A. Further, the WRKY1 and DEAD-box-RH27 were found to be associated with seed weight. Furthermore, the associations of PIF3 and pentatricopeptide-repeat-containing-gene with maturity and seed weight, and aldo–keto-reductase with flowering and maturity were revealed.ConclusionThis study offers insights into the genetic basis of key agronomic traits in rice bean, including flowering, maturity, and seed weight. The identified markers and associated candidate genes provide valuable resources for future exploration and targeted breeding, aiming to enhance the agronomic performance of rice bean cultivars. Notably, this research represents the first transcriptome-wide association study in pulse crop, uncovering the candidate genes for agronomically useful traits.

Identification and characterization of whole blood gene expression and splicing quantitative trait loci during early to mid-lactation of dairy cattle

BackgroundCharacterization of regulatory variants (e.g., gene expression quantitative trait loci, eQTL; gene splicing QTL, sQTL) is crucial for biologically interpreting molecular mechanisms underlying loci associated with complex traits. However, regulatory variants in dairy cattle, particularly in specific biological contexts (e.g., distinct lactation stages), remain largely unknown. In this study, we explored regulatory variants in whole blood samples collected during early to mid-lactation (22–150 days after calving) of 101 Holstein cows and analyzed them to decipher the regulatory mechanisms underlying complex traits in dairy cattle.ResultsWe identified 14,303 genes and 227,705 intron clusters expressed in the white blood cells of 101 cattle. The average heritability of gene expression and intron excision ratio explained by cis-SNPs is 0.28 ± 0.13 and 0.25 ± 0.13, respectively. We identified 23,485 SNP-gene expression pairs and 18,166 SNP-intron cluster pairs in dairy cattle during early to mid-lactation. Compared with the 2,380,457 cis-eQTLs reported to be present in blood in the Cattle Genotype-Tissue Expression atlas (CattleGTEx), only 6,114 cis-eQTLs (P < 0.05) were detected in the present study. By conducting colocalization analysis between cis-e/sQTL and the results of genome-wide association studies (GWAS) from four traits, we identified a cis-e/sQTL (rs109421300) of the DGAT1 gene that might be a key marker in early to mid-lactation for milk yield, fat yield, protein yield, and somatic cell score (PP4 > 0.6). Finally, transcriptome-wide association studies (TWAS) revealed certain genes (e.g., FAM83H and TBC1D17) whose expression in white blood cells was significantly (P < 0.05) associated with complex traits.ConclusionsThis study investigated the genetic regulation of gene expression and alternative splicing in dairy cows during early to mid-lactation and provided new insights into the regulatory mechanisms underlying complex traits of economic importance.

Abstract PO5-09-03: Multi-tissue transcriptome-wide association studies identified genes for intrinsic subtypes of breast cancer

Abstract Background: Though several genome-wide association studies (GWAS) of breast cancer (BC) have identified common variants which differ between intrinsic subtypes, the genes through which these variants act through to impact BC risk have not been fully established. Furthermore, transcriptome-wide association studies (TWAS) have thus far been primarily employed to identify genes associated with overall BC risk, while overlooking how the influence of common variation on gene expression may contribute to subtype-specific differences. Methods: In this study, we performed two complementary multi-tissue TWASs for each of the following BC intrinsic subtypes: Luminal A-like, Luminal B-like, Luminal B/HER2-negative-like, HER2-enriched-like, Triple Negative BC. These two approaches included 1) an expression-based approach that collated TWAS signals from expression quantitative trait loci (eQTLs) across multiple tissues using the Aggregated Cauchy Association Test (ACAT) and 2) a splicing-based approach that utilized two applications of ACAT to collate splicing QTLs (sQTLs) for a given gene in a tissue and then across tissues. To perform these two TWASs, we utilized e/sQTL models trained in 11 tissues from the Genotype-Tissue Expression Project including breast, ovary, uterus, vagina, EBV-transformed lymphocytes, whole blood, spleen, liver, subcutaneous adipose, visceral adipose, and cell-cultured fibroblasts. GWAS summary statistics were previously generated from 133,384 BC cases and 113,789 controls who were participants in the Breast Cancer Association Consortium (BCAC). We further performed our TWAS while conditioning e/sQTL effect sizes on nearby GWAS index SNPs. Additionally, we utilized gene-based fine-mapping of eQTLs and sQTL to identify candidate causal genes for each intrinsic subtype. Results: Overall, we identified 164 genes in 69 loci that were associated with Luminal A-like, 19 genes in 9 loci with Luminal B-like, 18 genes in 11 loci with Luminal B/HER2-negative-like, 10 genes in 7 loci with HER2-enriched-like, and 29 genes in 12 loci with TNBC. Among these genes, 17 genes had not been reported in previous TWAS of BC, and 140 genes, 1 gene, 2 genes, 2 genes, and 16 genes were uniquely associated with each of the intrinsic subtypes, respectively. Additionally, we identified one gene associated with Luminal A-like and one gene with Luminal B-like BC that were each in a locus located at least 1.4 Mb from published GWAS hits. Furthermore, we identified 106, 11, 10, 5, and 21 candidate causal genes for each of the intrinsic subtypes, respectively, that had a posterior inclusion probability of 0.9 in at least one e/sQTL. Conclusion: In summary, our multi-tissue TWAS corroborated previous GWAS loci for overall BC risk and intrinsic subtypes, while underscoring how common variation impacts BC etiology by modulating the expression and splicing of genes in multiple tissue types. Citation Format: James Li, Julian McClellan, Haoyu Zhang, Guimin Gao, Dezheng Huo. Multi-tissue transcriptome-wide association studies identified genes for intrinsic subtypes of breast cancer [abstract]. In: Proceedings of the 2023 San Antonio Breast Cancer Symposium; 2023 Dec 5-9; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2024;84(9 Suppl):Abstract nr PO5-09-03.

Using genome and transcriptome data from African-ancestry female participants to identify putative breast cancer susceptibility genes

African-ancestry (AA) participants are underrepresented in genetics research. Here, we conducted a transcriptome-wide association study (TWAS) in AA female participants to identify putative breast cancer susceptibility genes. We built genetic models to predict levels of gene expression, exon junction, and 3′ UTR alternative polyadenylation using genomic and transcriptomic data generated in normal breast tissues from 150 AA participants and then used these models to perform association analyses using genomic data from 18,034 cases and 22,104 controls. At Bonferroni-corrected P < 0.05, we identified six genes associated with breast cancer risk, including four genes not previously reported (CTD-3080P12.3, EN1, LINC01956 and NUP210L). Most of these genes showed a stronger association with risk of estrogen-receptor (ER) negative or triple-negative than ER-positive breast cancer. We also replicated the associations with 29 genes reported in previous TWAS at P < 0.05 (one-sided), providing further support for an association of these genes with breast cancer risk. Our study sheds new light on the genetic basis of breast cancer and highlights the value of conducting research in AA populations.

Integrative multi-omics analyses to identify the genetic and functional mechanisms underlying ovarian cancer risk regions

To identify credible causal risk variants (CCVs) associated with different histotypes of epithelial ovarian cancer (EOC), we performed genome-wide association analysis for 470,825 genotyped and 10,163,797 imputed SNPs in 25,981 EOC cases and 105,724 controls of European origin. We identified five histotype-specific EOC risk regions (p value <5×10-8) and confirmed previously reported associations for 27 risk regions. Conditional analyses identified an additional 11 signals independent of the primary signal at six risk regions (pvalue <10-5). Fine mapping identified 4,008 CCVs in these regions, of which 1,452 CCVs were located in ovarian cancer-related chromatin marks with significant enrichment in active enhancers, active promoters, and active regions for CCVs from each EOC histotype. Transcriptome-wide association and colocalization analyses across histotypes using tissue-specific and cross-tissue datasets identified 86 candidate susceptibility genes in known EOC risk regions and 32 genes in 23 additional genomic regions that may represent novel EOC risk loci (false discovery rate <0.05). Finally, by integrating genome-wide HiChIP interactome analysis with transcriptome-wide association study (TWAS), variant effect predictor, transcription factor ChIP-seq, and motifbreakR data, we identified candidate gene-CCV interactions at each locus. This included risk loci where TWAS identified one or more candidate susceptibility genes (e.g., HOXD-AS2, HOXD8, and HOXD3 at 2q31) and other loci where no candidate gene was identified (e.g., MYC and PVT1 at 8q24) by TWAS. In summary, this study describes a functional framework and provides a greater understanding of the biological significance of risk alleles and candidate gene targets at EOC susceptibility loci identified by a genome-wide association study.

Causal association and shared genetics between telomere length and COVID-19 outcomes: new evidence from the latest large-scale summary statistics

BackgroundObservational studies suggested that leukocyte telomere length (LTL) is shortened in COVID-19 patients. However, the genetic association and causality remained unknown. MethodsBased on the genome-wide association of LTL (N = 472,174) and COVID-19 phenotypes (N = 1,086,211-2,597,856), LDSC and SUPERGNOVA were used to estimate the genetic correlation. Cross-trait GWAS meta-analysis, colocalization, fine-mapping analysis and transcriptome-wide association study were conducted to explore the shared genetic etiology. Mendelian randomization (MR) was utilized to infer the causality. Upstream and downstream two-step MR was performed to investigate the potential mediating effects. ResultsLDSC identified a significant genetic association between LTL and all COVID-19 phenotypes (rG < 0, p < 0.05). Six significant regions were observed for LTL and COVID-19 susceptibility and hospitalization, respectively. Colocalization analysis found rs144204502, rs34517439, and rs56255908 were shared causal variants between LTL and COVID-19 phenotypes. Numerous biological pathways associated with LTL and COVID-19 outcomes were identified, mainly involved in immune related pathways. MR showed that longer LTL was significantly associated with a lower risk of COVID-19 severity (OR [95% CI] = 0.81 [0.71-0.92], p = 1.24×10-3) and suggestively associated with lower risks of COVID-19 susceptibility (OR [95% CI] = 0.96 [0.92-1.00], p = 3.44×10-2) and COVID-19 hospitalization (OR [95% CI] = 0.89 [0.80-0.98], p = 1.89×10-2). LTL partially mediated the effects of BMI, smoking, and education on COVID-19 outcomes. Furthermore, six proteins partially mediated the causality of LTL on COVID-19 outcomes, including BNDF, QPCT, FAS, MPO, SFTPB, and APOF. ConclusionsOur findings suggested that shorter LTL was genetically associated with higher risk of COVID-19 phenotypes, with shared genetic etiology and potential causality.

A multi-tissue, splicing-based joint transcriptome-wide association study identifies susceptibility genes for breast cancer.

Splicing-based transcriptome-wide association studies (splicing-TWASs) of breast cancer have the potential to identify susceptibility genes. However, existing splicing-TWASs test the association of individual excised introns in breast tissue only and thus have limited power to detect susceptibility genes. In this study, we performed a multi-tissue joint splicing-TWAS that integrated splicing-TWAS signals of multiple excised introns in each gene across 11 tissues that are potentially relevant to breast cancer risk. We utilized summary statistics from a meta-analysis that combined genome-wide association study (GWAS) results of 424,650 women of European ancestry. Splicing-level prediction models were trained in GTEx (v.8) data. We identified 240 genes by the multi-tissue joint splicing-TWAS at the Bonferroni-corrected significance level; in the tissue-specific splicing-TWAS that combined TWAS signals of excised introns in genes in breast tissue only, we identified nine additional significant genes. Of these 249 genes, 88 genes in 62 loci have not been reported by previous TWASs, and 17 genes in seven loci are at least 1 Mb away from published GWAS index variants. By comparing the results of our splicing-TWASs with previous gene-expression-based TWASs that used the same summary statistics and expression prediction models trained in the same reference panel, we found that 110 genes in 70 loci that are identified only by the splicing-TWASs. Our results showed that for many genes, expression quantitative trait loci (eQTL) did not show a significant impact on breast cancer risk, whereas splicing quantitative trait loci (sQTL) showed a strong impact through intron excision events.

Cis- and trans-eQTL TWASs of breast and ovarian cancer identify more than 100 susceptibility genes in the BCAC and OCAC consortia.

Transcriptome-wide association studies (TWASs) have investigated the role of genetically regulated transcriptional activity in the etiologies of breast and ovarian cancer. However, methods performed to date have focused on the regulatory effects of risk-associated SNPs thought to act in cis on a nearby target gene. With growing evidence for distal (trans) regulatory effects of variants on gene expression, we performed TWASs of breast and ovarian cancer using a Bayesian genome-wide TWAS method (BGW-TWAS) that considers effects of both cis- and trans-expression quantitative trait loci (eQTLs). We applied BGW-TWAS to whole-genome and RNA sequencing data in breast and ovarian tissues from the Genotype-Tissue Expression project to train expression imputation models. We applied these models to large-scale GWAS summary statistic data from the Breast Cancer and Ovarian Cancer Association Consortia to identify genes associated with risk of overall breast cancer, non-mucinous epithelial ovarian cancer, and 10 cancer subtypes. We identified 101 genes significantly associated with risk with breast cancer phenotypes and 8 with ovarian phenotypes. These loci include established risk genes and several novel candidate risk loci, such as ACAP3, whose associations are predominantly driven by trans-eQTLs. We replicated several associations using summary statistics from an independent GWAS of these cancer phenotypes. We further used genotype and expression data in normal and tumor breast tissue from the Cancer Genome Atlas to examine the performance of our trained expression imputation models. This work represents an in-depth look into the role of trans eQTLs in the complex molecular mechanisms underlying these diseases.

​Comprehensive mendelian randomization analysis of plasma proteomics to identify new therapeutic targets for the treatment of coronary heart disease and myocardial infarction.

Ischemic heart disease is one of the leading causes of mortality worldwide, and thus calls for development of more effective therapeutic strategies. This study aimed to identify potential therapeutic targets for coronary heart disease (CHD) and myocardial infarction (MI) by investigating the causal relationship between plasma proteins and these conditions. A two-sample Mendelian randomization (MR) study was performed to evaluate more than 1600 plasma proteins for their causal associations with CHD and MI. The MR findings were further confirmed through Bayesian colocalization, Summary-data-based Mendelian Randomization (SMR), and Transcriptome-Wide Association Studies (TWAS) analyses. Further analyses, including enrichment analysis, single-cell analysis, MR analysis of cardiovascular risk factors, phenome-wide Mendelian Randomization (Phe-MR), and protein-protein interaction (PPI) network construction were conducted to verify the roles of selected causal proteins. Thirteen proteins were causally associated with CHD, seven of which were also causal for MI. Among them, FES and PCSK9 were causal proteins for both diseases as determined by several analytical methods. PCSK9 was a risk factor of CHD (OR = 1.25, 95% CI: 1.13-1.38, P = 7.47E-06) and MI (OR = 1.36, 95% CI: 1.21-1.54, P = 2.30E-07), whereas FES was protective against CHD (OR = 0.68, 95% CI: 0.59-0.79, P = 6.40E-07) and MI (OR = 0.65, 95% CI: 0.54-0.77, P = 5.38E-07). Further validation through enrichment and single-cell analysis confirmed the causal effects of these proteins. Moreover, MR analysis of cardiovascular risk factors, Phe-MR, and PPI network provided insights into the potential drug development based on the proteins. This study investigated the causal pathways associated with CHD and MI, highlighting the protective and risk roles of FES and PCSK9, respectively. FES. Specifically, the results showed that these proteins are promising therapeutic targets for future drug development.

Integration of multi-omics summary data reveals the role of N6-methyladenosine in neuropsychiatric disorders.

N6-methyladenosine (m6A) methylation regulates gene expression/protein by influencing numerous aspects of mRNA metabolism and contributes to neuropsychiatric diseases. Here, we integrated multi-omics data and genome-wide association study summary data of schizophrenia (SCZ), bipolar disorder (BP), attention deficit hyperactivity disorder (ADHD), autism spectrum disorder (ASD), major depressive disorder (MDD), Alzheimer's disease (AD), and Parkinson's disease (PD) to reveal the role of m6A in neuropsychiatric disorders by using transcriptome-wide association study (TWAS) tool and Summary-data-based Mendelian randomization (SMR). Our investigation identified 86 m6A sites associated with seven neuropsychiatric diseases and then revealed 7881 associations between m6A sites and gene expressions. Based on these results, we discovered 916 significant m6A-gene associations involving 82 disease-related m6A sites and 606 genes. Further integrating the 58 disease-related genes from TWAS and SMR analysis, we obtained 61, 8, 7, 3, and 2 associations linking m6A-disease, m6A-gene, and gene-disease for SCZ, BP, AD, MDD, and PD separately. Functional analysis showed the m6A mapped genes were enriched in "response to stimulus" pathway. In addition, we also analyzed the effect of gene expression on m6A and the post-transcription effect of m6A on protein. Our study provided new insights into the genetic component of m6A in neuropsychiatric disorders and unveiled potential pathogenic mechanisms where m6A exerts influences on disease through gene expression/protein regulation.

Large-scale genome-wide association studies reveal the genetic causal etiology between air pollutants and autoimmune diseases

BackgroundEpidemiological evidence links a close correlation between long-term exposure to air pollutants and autoimmune diseases, while the causality remained unknown.MethodsTwo-sample Mendelian randomization (TSMR) was used to investigate the role of PM10, PM2.5, NO2, and NOX (N = 423,796–456,380) in 15 autoimmune diseases (N = 14,890–314,995) using data from large European GWASs including UKB, FINNGEN, IMSGC, and IPSCSG. Multivariable Mendelian randomization (MVMR) was conducted to investigate the direct effect of each air pollutant and the mediating role of common factors, including body mass index (BMI), alcohol consumption, smoking status, and household income. Transcriptome-wide association studies (TWAS), two-step MR, and colocalization analyses were performed to explore underlying mechanisms between air pollution and autoimmune diseases.ResultsIn TSMR, after correction of multiple testing, hypothyroidism was causally associated with higher exposure to NO2 [odds ratio (OR): 1.37, p = 9.08 × 10–4] and NOX [OR: 1.34, p = 2.86 × 10–3], ulcerative colitis (UC) was causally associated with higher exposure to NOX [OR: 2.24, p = 1.23 × 10–2] and PM2.5 [OR: 2.60, p = 5.96 × 10–3], rheumatoid arthritis was causally associated with higher exposure to NOX [OR: 1.72, p = 1.50 × 10–2], systemic lupus erythematosus was causally associated with higher exposure to NOX [OR: 4.92, p = 6.89 × 10–3], celiac disease was causally associated with lower exposure to NOX [OR: 0.14, p = 6.74 × 10–4] and PM2.5 [OR: 0.17, p = 3.18 × 10–3]. The risky effects of PM2.5 on UC remained significant in MVMR analyses after adjusting for other air pollutants. MVMR revealed several common mediators between air pollutants and autoimmune diseases. Transcriptional analysis identified specific gene transcripts and pathways interconnecting air pollutants and autoimmune diseases. Two-step MR revealed that POR, HSPA1B, and BRD2 might mediate from air pollutants to autoimmune diseases. POR pQTL (rs59882870, PPH4=1.00) strongly colocalized with autoimmune diseases.ConclusionThis research underscores the necessity of rigorous air pollutant surveillance within public health studies to curb the prevalence of autoimmune diseases.Graphical abstract (Built by the Biorender)

Large-scale cross-ancestry genome-wide meta-analysis of serum urate

Hyperuricemia is an essential causal risk factor for gout and is associated with cardiometabolic diseases. Given the limited contribution of East Asian ancestry to genome-wide association studies of serum urate, the genetic architecture of serum urate requires exploration. A large-scale cross-ancestry genome-wide association meta-analysis of 1,029,323 individuals and ancestry-specific meta-analysis identifies a total of 351 loci, including 17 previously unreported loci. The genetic architecture of serum urate control is similar between European and East Asian populations. A transcriptome-wide association study, enrichment analysis, and colocalization analysis in relevant tissues identify candidate serum urate-associated genes, including CTBP1, SKIV2L, and WWP2. A phenome-wide association study using polygenic risk scores identifies serum urate-correlated diseases including heart failure and hypertension. Mendelian randomization and mediation analyses show that serum urate-associated genes might have a causal relationship with serum urate-correlated diseases via mediation effects. This study elucidates our understanding of the genetic architecture of serum urate control.

Prioritization of therapeutic targets for cancers using integrative multi-omics analysis

BackgroundThe integration of transcriptomic, proteomic, druggable genetic and metabolomic association studies facilitated a comprehensive investigation of molecular features and shared pathways for cancers’ development and progression.MethodsComprehensive approaches consisting of transcriptome-wide association studies (TWAS), proteome-wide association studies (PWAS), summary-data-based Mendelian randomization (SMR) and MR were performed to identify genes significantly associated with cancers. The results identified in above analyzes were subsequently involved in phenotype scanning and enrichment analyzes to explore the possible health effects and shared pathways. Additionally, we also conducted MR analysis   to investigate metabolic pathways related to cancers.ResultsTotally 24 genes (18 transcriptomic, 1 proteomic and 5 druggable genetic) showed significant associations with cancers risk. All genes identified in multiple methods were mainly enriched in nuclear factor erythroid 2-related factor 2 (NRF2) pathway. Additionally, biosynthesis of ubiquinol and urate were found to play an important role in gastrointestinal tumors.ConclusionsA set of putatively causal genes and pathways relevant to cancers were identified in this study, shedding light on the shared biological processes for tumorigenesis and providing compelling genetic evidence to prioritize anti-cancer drugs development.

Multi-tissue transcriptome-wide association study identifies novel candidate susceptibility genes for cataract

IntroductionCataract is the leading cause of blindness among the elderly worldwide. Twin and family studies support an important role for genetic factors in cataract susceptibility with heritability estimates up to 58%. To date, 55 loci for cataract have been identified by genome-wide association studies (GWAS), however, much work remains to identify the causal genes. Here, we conducted a transcriptome-wide association study (TWAS) of cataract to prioritize causal genes and identify novel ones, and examine the impact of their expression.MethodsWe performed tissue-specific and multi-tissue TWAS analyses to assess associations between imputed gene expression from 54 tissues (including 49 from the Genotype Tissue Expression (GTEx) Project v8) with cataract using FUSION software. Meta-analyzed GWAS summary statistics from 59,944 cataract cases and 478,571 controls, all of European ancestry and from two cohorts (GERA and UK Biobank) were used. We then examined the expression of the novel genes in the lens tissue using the iSyTE database.ResultsAcross tissue-specific and multi-tissue analyses, we identified 99 genes for which genetically predicted gene expression was associated with cataract after correcting for multiple testing. Of these 99 genes, 20 (AC007773.1, ANKH, ASIP, ATP13A2, CAPZB, CEP95, COQ6, CREB1, CROCC, DDX5, EFEMP1, EIF2S2, ESRRB, GOSR2, HERC4, INSRR, NIPSNAP2, PICALM, SENP3, and SH3YL1) did not overlap with previously reported cataract-associated loci. Tissue-specific analysis identified 202 significant gene-tissue associations for cataract, of which 166 (82.2%), representing 9 unique genes, were attributed to the previously reported 11q13.3 locus. Tissue-enrichment analysis revealed that gastrointestinal tissues represented one of the highest proportions of the Bonferroni-significant gene-tissue associations (21.3%). Moreover, this gastrointestinal tissue type was the only anatomical category significantly enriched in our results, after correcting for the number of tissue donors and imputable genes for each reference panel. Finally, most of the novel cataract genes (e.g., Capzb) were robustly expressed in iSyTE lens data.DiscussionOur results provide evidence of the utility of imputation-based TWAS approaches to characterize known GWAS risk loci and identify novel candidate genes that may increase our understanding of cataract etiology. Our findings also highlight the fact that expression of genes associated with cataract susceptibility is not necessarily restricted to lens tissue.

Investigating the role of common cis-regulatory variants in modifying penetrance of putatively damaging, inherited variants in severe neurodevelopmental disorders

Recent work has revealed an important role for rare, incompletely penetrant inherited coding variants in neurodevelopmental disorders (NDDs). Additionally, we have previously shown that common variants contribute to risk for rare NDDs. Here, we investigate whether common variants exert their effects by modifying gene expression, using multi-cis-expression quantitative trait loci (cis-eQTL) prediction models. We first performed a transcriptome-wide association study for NDDs using 6987 probands from the Deciphering Developmental Disorders (DDD) study and 9720 controls, and found one gene, RAB2A, that passed multiple testing correction (p = 6.7 × 10–7). We then investigated whether cis-eQTLs modify the penetrance of putatively damaging, rare coding variants inherited by NDD probands from their unaffected parents in a set of 1700 trios. We found no evidence that unaffected parents transmitting putatively damaging coding variants had higher genetically-predicted expression of the variant-harboring gene than their child. In probands carrying putatively damaging variants in constrained genes, the genetically-predicted expression of these genes in blood was lower than in controls (p = 2.7 × 10–3). However, results for proband-control comparisons were inconsistent across different sets of genes, variant filters and tissues. We find limited evidence that common cis-eQTLs modify penetrance of rare coding variants in a large cohort of NDD probands.

Shared Genetic Determinants of Schizophrenia and Autism Spectrum Disorder Implicate Opposite Risk Patterns: A Genome-Wide Analysis of Common Variants.

The synaptic pruning hypothesis posits that schizophrenia (SCZ) and autism spectrum disorder (ASD) may represent opposite ends of neurodevelopmental disorders: individuals with ASD exhibit an overabundance of synapses and connections while SCZ was characterized by excessive pruning of synapses and a reduction. Given the strong genetic predisposition of both disorders, we propose a shared genetic component, with certain loci having differential regulatory impacts. Genome-Wide single nucleotide polymorphism (SNP) data of European descent from SCZ (N cases = 53386, N controls = 77258) and ASD (N cases = 18381, N controls = 27969) were analyzed. We used genetic correlation, bivariate causal mixture model, conditional false discovery rate method, colocalization, Transcriptome-Wide Association Study (TWAS), and Phenome-Wide Association Study (PheWAS) to investigate the genetic overlap and gene expression pattern. We found a positive genetic correlation between SCZ and ASD (rg = .26, SE = 0.01, P = 7.87e-14), with 11 genomic loci jointly influencing both conditions (conjFDR <0.05). Functional analysis highlights a significant enrichment of shared genes during early to mid-fetal developmental stages. A notable genetic region on chromosome 17q21.31 (lead SNP rs2696609) showed strong evidence of colocalization (PP.H4.abf = 0.85). This SNP rs2696609 is linked to many imaging-derived brain phenotypes. TWAS indicated opposing gene expression patterns (primarily pseudogenes and long noncoding RNAs [lncRNAs]) for ASD and SCZ in the 17q21.31 region and some genes (LRRC37A4P, LINC02210, and DND1P1) exhibit considerable variation in the cerebellum across the lifespan. Our findings support a shared genetic basis for SCZ and ASD. A common genetic variant, rs2696609, located in the Chr17q21.31 locus, may exert differential risk regulation on SCZ and ASD by altering brain structure. Future studies should focus on the role of pseudogenes, lncRNAs, and cerebellum in synaptic pruning and neurodevelopmental disorders.

A new test for trait mean and variance detects unreported loci for blood-pressure variation

Variability in quantitative traits has clinical, ecological, and evolutionary significance. Most genetic variants identified for complex quantitative traits have only a detectable effect on the mean of trait. We have developed the mean-variance test (MVtest) to simultaneously model the mean and log-variance of a quantitative trait as functions of genotypes and covariates by using estimating equations. The advantages of MVtest include the facts that it can detect effect modification, that multiple testing can follow conventional thresholds, that it is robust to non-normal outcomes, and that association statistics can be meta-analyzed. In simulations, we show control of type I error of MVtest over several alternatives. We identified 51 and 37 previously unreported associations for effects on blood-pressure variance and mean, respectively, in the UK Biobank. Transcriptome-wide association studies revealed 633 significant unique gene associations with blood-pressure mean variance. MVtest is broadly applicable to studies of complex quantitative traits and provides an important opportunity to detect novel loci.

Identification of candidate genes associated with milk production and mastitis based on transcriptome-wide association study.

Genetic research for the assessment of mastitis and milk production traits simultaneously has a long history. The main issue that arises in this context is the known existence of a positive correlation between the risk of mastitis and lactation performance due to selection. The transcriptome-wide association study (TWAS) approach endeavors to combine the expression quantitative trait loci and genome-wide association study summary statistics to decode complex traits or diseases. Accordingly, we used the farmgtex project results as a complete bovine database for mastitis and milk production. The results of colocalization and TWAS approaches were used for the detection of functional associated candidate genes with milk production and mastitis traits on multiple tissue-based transcriptome records. Also, we used the david database for gene ontology to identify significant terms and associated genes. For the identification of interaction networks, the genemania and string databases were used. Also, the available z-scores in TWAS results were used for the calculation of the correlation between tissues. Therefore, the present results confirm that LYNX1, DGAT1, C14H8orf33, and LY6E were identified as significant genes associated with milk production in eight, six, five, and five tissues, respectively. Also, FBXL6 was detected as a significant gene associated with mastitis trait. CLN3 and ZNF34 genes emerged via both the colocalization and TWAS approaches as significant genes for milk production trait. It is expected that TWAS and colocalization can improve our perception of the potential health status control mechanism in high-yielding dairy cows.

The broad impact of cell death genes on the human disease phenome

Cell death mediated by genetically defined signaling pathways influences the health and dynamics of all tissues, however the tissue specificity of cell death pathways and the relationships between these pathways and human disease are not well understood. We analyzed the expression profiles of an array of 44 cell death genes involved in apoptosis, necroptosis, and pyroptosis cell death pathways across 49 human tissues from GTEx, to elucidate the landscape of cell death gene expression across human tissues, and the relationship between tissue-specific genetically determined expression and the human phenome. We uncovered unique cell death gene expression profiles across tissue types, suggesting there are physiologically distinct cell death programs in different tissues. Using summary statistics-based transcriptome wide association studies (TWAS) on human traits in the UK Biobank (n ~ 500,000), we evaluated 513 traits encompassing ICD-10 defined diagnoses and laboratory-derived traits. Our analysis revealed hundreds of significant (FDR < 0.05) associations between genetically regulated cell death gene expression and an array of human phenotypes encompassing both clinical diagnoses and hematologic parameters, which were independently validated in another large-scale DNA biobank (BioVU) at Vanderbilt University Medical Center (n = 94,474) with matching phenotypes. Cell death genes were highly enriched for significant associations with blood traits versus non-cell-death genes, with apoptosis-associated genes enriched for leukocyte and platelet traits. Our findings are also concordant with independently published studies (e.g. associations between BCL2L11/BIM expression and platelet & lymphocyte counts). Overall, these results suggest that cell death genes play distinct roles in their contribution to human phenotypes, and that cell death genes influence a diverse array of human traits.